Method and Apparatus for the Stabilization of Spectrometric Transducers

ABSTRACT

The invention provides a method for stabilising a spectrometric transducer for an optical spectrum measuring instrument by obtaining an instantaneous central wavelength of a thermally controlled tunable filter of the instrument, calibrating for a selected ambient temperature, determining a heat variance of the filter and controlling the filter to compensate for heat variance

This application is a continuation or U.S. patent application Ser. No. 12/901,238 filed Oct. 8, 2011 which is a continuation of U.S. patent application Ser. No. 12/318,040 filed on Dec. 19, 2008 which is a continuation of U.S. patent application Ser. No. 11/754,781 filed on May 29, 2007 which is a continuation of U.S. patent application Ser. No. 11/038,490 on Jan. 21, 2005 which claims priority from U.S. Provisional Patent Application No. 60/537,919 filed on Jan. 22, 2004; the contents which are incorporated herein by reference.

FIELD OF THE INVENTION

The present invention relates to a method for the stabilization of spectrometric transducers.

BRIEF DESCRIPTION OF THE DRAWINGS

An embodiment of the present invention will now be described by way of example only with reference to the accompanying drawings, in which:

FIG. 1 is a schematic diagram of an optical spectrum measurement instrument.

FIG. 2 is a schematic representation of a portion of the apparatus shown in FIG. 1.

FIG. 3 is a flow diagram of a thermal stabilization method for a spectrometric transducer.

FIG. 4 is a plot showing the variation of controlling current for a particular set of parameters.

DETAILED DESCRIPTION OF THE INVENTION

In FIG. 1 there is shown an optical spectrum measurement instrument 20 for analyzing optical signal 12. The optical signal 12 is sampled by sampler 14, which provides optical spectrum x(λ) 16 to the optical spectrum measurement instrument 20. The optical spectrum measurement instrument 20 includes spectrometric transducer 30, including a tunable filter 32 and photodiode 34, and Digital Signal Processor (DSP) 40. Optical spectrum x(λ) 16 is separated into its optical components by tunable filter 32, after which its components are converted from optical to current by photodiodes 34, resulting in an electrical spectrometric data representation y(t) 36 of optical spectrum x(λ) 16. The spectrometric data representation y(t) 36 is then usedby the DSP 40 to compute an estimate {tilde over (x)}(λ) 42 of the optical signal 12 optical spectrum, which is then provided to user display 50. Optionally, the DSP 40 may also compute various parameters relating to the optical signal 12 optical spectrum such as power measurements, OSNR, BER, etc. or may implement signal analysis or reconstruction algorithms. The response of a filter to an input optical signal whose spectrum is x(λ) 16, i.e. the output current of the photodiode 34, may be modelled by an integral operator of the form:

$\begin{matrix} {{y(t)} = {F_{y}\left\lbrack {\int_{- \infty}^{+ \infty}{{g\left( {\lambda,t} \right)} \cdot {x(\lambda)} \cdot {\lambda}}} \right\rbrack}} & {{Equation}\mspace{14mu} 1} \end{matrix}$

where λ is wavelength, t is time, g(λ, t) is the filter response function, and F_(y)[•] is a slightly nonlinear scalar function of a scalar variable. For a monochromatic input signal, whose spectrum is x(λ)=δ(λ−l), this model responds with:

$\begin{matrix} {{y(t)} = {{F_{y}\left\lbrack {\int_{- \infty}^{+ \infty}{{g\left( {\lambda,t} \right)} \cdot {\delta \left( {\lambda - l} \right)} \cdot {\lambda}}} \right\rbrack} = {F_{y}\left\lbrack {g\left( {l,t} \right)} \right\rbrack}}} & {{Equation}\mspace{14mu} 2} \end{matrix}$

and to the input signal whose spectrum is flat, x(λ)=1, it responds with:

$\begin{matrix} {{y(t)} = {{F_{y}\left\lbrack {\int_{- \infty}^{+ \infty}{{g\left( {\lambda,t} \right)}\  \cdot {\lambda}}} \right\rbrack} \equiv {{F_{y}\left\lbrack {g_{0}(t)} \right\rbrack}.}}} & {{Equation}\mspace{14mu} 3} \end{matrix}$

The response g(l, t), defined by Equation 2, is a gaussoid-like function with a maximum changing monotonically with the wavelength l. The function of time g₀(t), defined by Equatiou 3, characterizes the amplitude variability of the filter response along the wavelength axis.

As shown in Figure 2, the instantaneous central wavelength of the filter is externally controlled by the current i(t) 38. In the preferred embodiment, the filter 32 is a thermally controlled filter in which the central wavelength depends on its internal temperature. The current i(t) 38 drives a variable resistance heater incorporated in the filter. Such filters utilize the characteristic that as the temperature on tunable filters varies, such as a particular optical material's change in its index of refraction n with temperature T, so its center wavelength varies, resulting in thermally-controlled tunable filters. The temperature of the filter depends on the heat supplied from the resistence heater and the heat loss/gain to/from the ambient environment. To compensate for the heat loss/gain, a signal containing information about the ambient temperature is used to compensate the current i(t) 38 supplied to the heater. The current is adjusted by the controller that implements a control algorithm to maintain the filter at the desired center frequency.

Assume that the function λ=F₀(u) is the result ofwayelength calibration for a selected ambient temperature T₀ represented by the voltage U₀. Assume, moreover, that the input signal is generated by a battery of lasers whose central wavelengths are uniformly distributed in the interval [λ_(min), λ_(max)]:

$\begin{matrix} {{\lambda_{j} = {{\lambda_{\min} + {\frac{\left( {j - 1} \right)}{J - 1}\left( {\lambda_{\max} - \lambda_{\min}} \right)\mspace{14mu} {for}\mspace{14mu} j}} = 1}},\ldots \mspace{14mu},{J.}} & {{Equation}\mspace{14mu} 4} \end{matrix}$

Then, the maxima y_(0,j) ^(max) of the response of the stabilized filter to such a signal may be easily identified together with the corresponding values of the voltage u(t): u_(0,j)=F₀ ⁻¹(λ_(j)) for j=1, . . . , J. The practical purpose of controlling the current i(t) 38 may be now formulated as follows. For an arbitrary ambient temperature T, find i(t) such that the maxima of the stabilized filter response remain as close as possible to the maxima determined for T₀, i.e. the coordinates of the maxima for the ambient temperature T, the values u₁, . . . , u_(j) of the voltage u(t) and the values y₁ ^(max), . . . , y_(j) ^(max) of the output signal y(t), satisfy the condition:

w[(u ₁ −u _(0,1))²+. . . , +(u _(j) −u _(0,j))²]+(1−w)[(y ₁ ^(max) −y _(0,1) ^(max))²+. . . +(y _(j) ^(max) −y _(0,j) ^(max))²]→MIN   Equation 5

where w∈[0, 1] is a weighing factor.

The task of control may be significantly simplified (at the price of the sub-optimality of the solution) by an appropriate parameterization of the heating current i(t). An example of such a parameterization is defined by the formula:

$\begin{matrix} {{i(t)} = \left\{ \begin{matrix} 0 & {{{for}\mspace{14mu} t} < 0} \\ i_{1} & {{{for}\mspace{14mu} t} \in \left\lbrack {0,t_{2}} \right\rbrack} \\ {{\frac{i_{3} - i_{2}}{t_{3} - t_{2}}\left( {t - t_{2}} \right)} + i_{2}} & {{{for}\mspace{14mu} t} \in \left\lbrack {t_{2},t_{3}} \right\rbrack} \end{matrix} \right.} & {{Equation}\mspace{14mu} 6} \end{matrix}$

where i₁, i₂, i₃, t₂, and t₃ are parameters of the current to be optimized by minimization of the left-hand side of Equation 5. The parameterization of controlled current enables one to use an empirical procedure of optimization that does not require any explicit reference to the mathematical model of the stabilized filter. This procedure is implemented in the controller and is depicted by the flow chart shown in FIG 3. The sequence of steps composing the procedure is indicated by the sequence of blocks 62 to 70.

In block 62, for the selected ambient temperature T₀, the values of i₁, i₂, i₃, t₂ and t₃ are chosen so as to produce a relatively uniform distribution of u_(0,j)=F₀ ⁻¹(λ_(j)) for j=1, . . . , J.

Then in block 64, on the basis of measurements performed for the same ambient temperature, the matrix S_(y) of the sensitivity of the maxima y₁ ^(max), . . . , y_(j) ^(max) is computed for a small change in the ambient temperature ΔT and for small changes Δi₁, Δi₂, Δi₃, Δt₂ and Δt₃ in the parameters i₁, i₂, i₃, t₂ and t₃.

Similarly, in block 66, on the basis of measurements performed for the same ambient temperature, the matrix S_(u) of the sensitivity of the corresponding maxima voltage values u₁, . . . , u_(j) is computed for a small change in the ambient temperature ΔT and for small changes Δi₁, Δi₂, Δi₃, Δt₂ and Δt₃ in the parameters i₁, i₂, i₃, t₂ and t₃.

In block 68, the following minimization problem is solved:

$\begin{matrix} \left. {{w{{S_{u}\begin{bmatrix} {\Delta \; T} \\ {\Delta \; i_{1}} \\ {\Delta \; i_{2}} \\ {\Delta \; i_{3}} \\ {\Delta \; t_{2}} \\ {\Delta \; t_{3}} \end{bmatrix}}}^{2}} + {\left( {1 - w} \right){{S_{y}\begin{bmatrix} {\Delta \; T} \\ {\Delta \; i_{1}} \\ {\Delta \; i_{2}} \\ {\Delta \; i_{3}} \\ {\Delta \; t_{2}} \\ {\Delta \; t_{3}} \end{bmatrix}}}^{2}}}\rightarrow{{MIN}.} \right. & {{Equation}\mspace{14mu} 7} \end{matrix}$

with respect to Δi₁, Δi₂, Δi₃, Δt₂ and Δt₃ for an assumed (sufficiently small) increment ΔT of the ambient temperature.

Following which, in block 70, the ambient temperature is changed to the value T=T₀+ΔT; the new values of i₁, i₂, i₃, t₂ and t₃ are computed using increments Δi₁, Δi₂, Δi₃, Δt₂ and Δt₃ determined in block 68; i₁, i₂, i₃, t₂ and t₃ and are empirically corrected in such a way as to satisfy the condition defined by Equation 5.

Finally, blocks 62 to 70 are repeated itcratively as to cover the whole range of ambient temperatures the stabilized filter is assumed to operate in.

It should be noted that the whole process of optimization is subject to the constraint concerning the admissible values of current and heating time.

A further refinement of the thermal stabilization of the stabilized filter is possible during the software pre-processing of the data provided by the stabilized filter using the above described hardware means. The residual instabilities may be characterized during general calibration of the stabilized filter, and the results of this characterization may next be used for correction of the raw data before their pre-processing.

A closer empirical study of the current i(t) control based on Equation 6 provides the following:

-   -   t₂ and t₃ may be fixed to the values 1 ms and 5 ms,         respectively;     -   abrupt changes of current i(t) should be avoided;     -   a convex control in the interval [t₂, t₃] would be desirable.

Consequently, a second example of the current i(t) parameterization may been designed.

It is defined by the formula:

$\begin{matrix} {{i(t)} = \left\{ \begin{matrix} 0 & {{{for}\mspace{14mu} t} < 0} \\ {{a_{1}t^{2}} + {b_{1}t} + c_{1}} & {{{for}\mspace{14mu} t} \in \left\lbrack {0,1} \right\rbrack} \\ {{a_{2}\left( {t - 1} \right)}^{2} + {b_{2}\left( {t - 1} \right)} + c_{2}} & {{{for}\mspace{14mu} t} \in \left\lbrack {1,5} \right\rbrack} \end{matrix} \right.} & {{Equation}\mspace{14mu} 8} \end{matrix}$

where the parameters a₁, b₁, c₁ and a₂, b₂, c₂ should satisfy the following conditions:

i(0)=I ₀ , i(1)=I ₁ , i(3)=I ₃, and i(5)=I ₅   Equation 9

The solution of the above algebraic problem has the form:

$\begin{matrix} {\begin{bmatrix} a_{1} \\ b_{1} \\ a_{2} \\ b_{2} \end{bmatrix} = {\begin{bmatrix} 1.000 & {- 1.750} & 1.000 & {- 0.250} \\ {- 2.000} & 2.750 & {- 1.000} & 0.250 \\ 0.000 & 0.125 & {- 0.250} & 0.125 \\ 0.000 & {- 0.750} & 1.000 & {- 0.250} \end{bmatrix} \cdot \begin{bmatrix} I_{0} \\ I_{1} \\ I_{3} \\ I_{5} \end{bmatrix}}} & {{Equation}\mspace{14mu} 10} \end{matrix}$

and c₁=I₀, c₂=I₁. An example of the current i(t) generated according to the above formula is shown in FIG. 4 for I₀=100 mA, I₁=50 mA, I₃=120 mA, and I₅=150 mA.

The components of a spectrometer or the device in which it is used such as an OPM are subject to aging. Consequently, the parameters of the OPM drift in time. In particular, the absolute accuracy of wavelength estimation is deteriorating. Taking into account that the contemporary DWDM transmitters contain highly stable lasers, one may use the time series of wavelength estimates as the basis for compensation of the time drift of the OPM.

The idea of using time series of wavelength estimates provided by the OPM for time drift compensation of this OPM is based on an assumption that the average deviation of the central wavelength of the laser signal used for this purpose is close to zero, i.e. there is no systematic evolution of this wavelength over time. Consequently, the average deviation of the central wavelength as computed by the OPM should be expected to also be close to zero. If not, this average deviation may be used to model the OPM's drift and provide a way of stabilizing the OPM by compensating for this time drift.

A wavelength such as described above may be modelled by means of a stochastic process:

l(t)={circumflex over (l)}+δl(t)  Equation 11

where t is time, {circumflex over (l)} is the central wavelength value according to the ITU grid, and δl(t) is a stochastic process modelling the wavelength deviation from the ITU-grid value. The latter process is assumed to be zero-mean and stationary. Consequently, the time sampling of l(t) at equidistant time points, t₁, . . . , t_(N), yields the vector of random variables:

l=[l(t₁) . . . l(t _(N))]^(T)   Equation 12

such that:

$\begin{matrix} {{{E\left\lbrack \underset{\_}{1} \right\rbrack} = {{\left\lbrack {\text{?}\; \ldots \mspace{14mu} \text{?}} \right\rbrack^{T}\mspace{14mu} {and}\mspace{14mu} {{Cov}\left\lbrack \underset{\_}{1} \right\rbrack}} = \begin{bmatrix} \sigma^{2} & \ldots & c_{1,N} \\ \vdots & \ddots & \vdots \\ c_{N,1} & \ldots & \sigma^{2} \end{bmatrix}^{T}}}{\text{?}\text{indicates text missing or illegible when filed}}} & {{Equation}\mspace{14mu} 13} \end{matrix}$

Under the above assumptions, a result of central wavelength measurement, provided by the OPM drifting in time, may be modelled with:

{circumflex over (l)}(t _(n))= l(t _(n))+Δl(t _(n))={circumflex over (l)}+Δl(t _(n))+δl(t _(n)) for n=1, . . . , N   Equation 14

where Δl(t_(n)) is the time drift of wavelength to be estimated on the basis of the realizations {circumflex over (l)}(t_(n)) of the random variables {circumflex over (l)}(t_(n)). The vector {circumflex over (I)}=[{circumflex over (l)}(t

). . . {circumflex over (l)}(t_(N))]^(T) has the following statistical properties:

$\begin{matrix} {{{E\lbrack i\rbrack} = \left\lbrack {i + {\Delta \; {i\left( t_{1} \right)}\mspace{14mu} \ldots \mspace{14mu} i} + {\Delta \; {i\left( t_{N} \right)}}} \right\rbrack^{T}}{and}{{{Cov}\left\lbrack \underset{\_}{i} \right\rbrack} = {\begin{bmatrix} \sigma^{2} & \ldots & c_{1,N} \\ \vdots & \ddots & \vdots \\ c_{1,N} & \ldots & \sigma^{2} \end{bmatrix}^{T} \equiv \Sigma}}} & {{Equation}\mspace{14mu} 15} \end{matrix}$

Assuming that the solution is to be approximated by a linear combination:

$\begin{matrix} {{{\Delta \text{?}(t)} = {\sum\limits_{j = 1}^{j}\; {P_{j}{\phi_{j}(t)}}}}{\text{?}\text{indicates text missing or illegible when filed}}} & {{Equation}\mspace{14mu} 16} \end{matrix}$

of known functions φ_(j)(t), such as, for example, a polynomial, an orthogonal polynomial, a trigonometric polynomial, a b-spline. etc, with unknown coefficients p_(j) forming the vector p=[p₁ . . . p_(j)]^(T). This vector is to be estimated on the oasis of the approximate equations:

$\begin{matrix} {\mspace{20mu} {{{\sum\limits_{j = 1}^{J}\; {p_{j}{\phi_{j}\left( t_{n} \right)}}} \equiv {{\text{?}\left( t_{n} \right)} - \text{?}}}\mspace{20mu} {{{{for}\mspace{14mu} n} = 1},\ldots \mspace{14mu},N}{\text{?}\text{indicates text missing or illegible when filed}}}} & {{Equation}\mspace{14mu} 17} \end{matrix}$

that may be given a matrix form:

Φ_(p)≅Δ{circumflex over (l)}  Equation 18

where:

$\begin{matrix} {\mspace{20mu} {{\Phi = \begin{bmatrix} {\phi_{1}\left( t_{1} \right)} & \ldots & {\phi_{j}\left( t_{1} \right)} \\ \vdots & \ddots & \vdots \\ {\phi_{1}\left( t_{N} \right)} & \ldots & {\phi_{j}\left( t_{N} \right)} \end{bmatrix}}\mspace{20mu} {and}\mspace{20mu} {{\Delta \text{?}} = {{\begin{bmatrix} \text{?} \\ \vdots \\ \text{?} \end{bmatrix}.\text{?}}\text{indicates text missing or illegible when filed}}}}} & {{Equation}\mspace{14mu} 19} \end{matrix}$

Under an assumption that there is no correlation between consecutive samples provided by the OPM, i.e. all c_(n,v)=0, the LS solution of this equation has the form:

{circumflex over (p)}=(Φ^(T)Φ)⁻¹Φ^(T)Δ{circumflex over (l)}  Equation 20

and the corresponding covariance matrix is:

Cov[{circumflex over (p)}]=(Φ^(T)Φ)⁻¹σ²   Equation 21

If the correlation between samples cannot be neglected, then the solution takes on the form:

{circumflex over (p)}=(Φ^(T)Σ⁻¹Φ)⁻¹Φ^(T)Σ⁻¹Δ{circumflex over (l)}  Equation 22

and the corresponding covariance matrix is:

Cov[{circumflex over (p)}]=(Φ^(T)Σ⁻¹Φ)⁻¹   Equation 23

If the samples may be considered uncorrected, then estimates {circumflex over (p)}₁, . . . , {circumflex over (p)}_(K) of the parameters p₁, . . . , p_(K), characterizing the drift of K DWDM channels, may be obtained using the LS method, described in the previous section, in an integrated numerical process:

[{circumflex over (p)}₁ . . . {circumflex over (p)}_(K)]=(Φ^(T)Φ)⁻¹Φ^(T)[Δ{circumflex over (l)}₁ . . . Δ{circumflex over (l)}_(K)]  Equation 24

where Δ{circumflex over (l)}_(k) is the vector of deviations of the results of measurements of the central wavelength in the kth channel from the ITU-grid value of this wavelengthth. In general, the channels may differ in the laser wavelength deviation; the corresponding variances σ₁ ², . . . , σ_(K) ² may be estimated

according to the formula:

$\begin{matrix} {{{\hat{\sigma}}_{k}^{2} = {\frac{1}{N - J}{{{\Phi {\hat{p}}_{k}} - {\Delta {\hat{l}}_{k}}}}_{2}^{2}}}{for}{{k = 1},\ldots \mspace{14mu},K}} & {{Equation}\mspace{14mu} 25} \end{matrix}$

Then the uncertainty-based weighing may be applied to obtain the solution:

$\begin{matrix} {\hat{p} = \frac{\sum\limits_{k = 1}^{K}\; {\frac{1}{\sigma_{k}^{2}}{\hat{p}}_{k}}}{\sum\limits_{k = 1}^{K}\; \frac{1}{\sigma_{k}^{2}}}} & {{Equation}\mspace{14mu} 26} \end{matrix}$

characterizing the wavelength-averaged drift of the OPM over time.

If the maximum deviation of the central wavelength of the laser signal used for the OPM drift compensation is negligible with respect to the assumed maximum error of wavelength measurements performed by the OPM, then the time drift of the OPM may be corrected on the basis of a single result of measurement.

As well, if the standard deviation of the central wavelength averaged over a time interval ΔT is negligible with respect to the assumed standard deviation of the error of wavelength measurements performed by the OPM, and the OPM drift during the time interval ΔT is negligible, then the time drift of OPM may be corrected on the basis of the average of results of measurement performed during ΔT.

Although the present invention has been described by way of a particular embodiment and examples thereof, it should be noted that it will be apparent to persons skilled in the art that modifications may be applied to the present particular embodiment without departing from the scope of the present invention. 

1. A method for stabilizing a spectrometric transducer for an optical spectrum measurement instrument comprising: obtaining an instantaneous central wavelength of a thermally controlled tunable filter of the instrument; calibrating for a selected ambient temperature: determining a heat variance of the filter; and controlling the filter to compensate for the heat variance. 